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ABSTRACT 


The general equations of motion for a rigid body are 
derived in cylindrical coordinates by Lagrangian dynamics 
and used to model the motion of an airplane in a steady spin. 
meeor Simplification, the equations are cast into a form 
MeLlizing conventional aerodynamic derivatives along with 
other derivatives which may be Significant in Spins. An 
iterative numerical solution procedure is outlined which 
should simplify the problem of solving the nonlinear dif- 
ferential equations, and relationships between the Euler 
memes used in the equations and the more familiar ordered 
Paamor pitch, roll, and yaw are derived to permit computer 
mieuc and output of orientation to be more easily visualized. 
The inverse problem is alSo considered, and equations are 
derived to translate spin test data into the six Lagrange 
coordinates used so that computer calculations may be checked 
. against test results and actual aerodynamic forces can be 
compared with computed values. The paper concludes with a 
discussion of the implications of this spin analysis along 


with applications and extensions. 
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TABLE OF SYMBOLS 


The definitions of the symbols used in this paper are as 


follows: 


Cry. Cry, Crs 


Cur, Cuy Cry 


J. , Gy, Jz 
I, 1, 4: 
M 

M, My Mz 


n 
re 4,7 


WM 


Center of gravity 
Mean aerodynamic chord, feet 


Aeredynamic Porce coefficient along the 
Subscripted axis, dimensionless 


Aerodynamic moment coefficient about the 
subscripted axis, dimensionless 


Aerodynamic force along the subscripted 
axis, lbs. 


Lagrange generalized non-conservative 
force, lbs. 


Lagrange generalized force in the subscripted 
coorcimares, Lb. or ft — lbs? 


Acceleration of gravity, ft/sec“ 


Airplane cg acceleration components along 
the subscripted axis, ft/sec2 


Airplane penerhey MOMENES sol eLNnew Caney, 
Slug - ft 


Airplane mass, slugs 
Moments about the subscripted axis, ft - lbs. 
Sequential subscript, integer 


Angular rates of rotation about the airplane 
X, Y, and Z principal axes, respectively 


Lagrange generalized coordinate 

Spin radius, feet 

Airplane wing area, ger 

Lagrange kinetic energy term, ft - lbs. 


Time, seconds 





At 


Meg 
Vann, Vauy, Vang 
x ¥2 

X,, Ys, 2 

Z. 

Z, 


7 


Increment of time, seconds 


Lagrange petential energy Germ, 7G) eos. 


Airplane cg velocity ft/sec 


Airplane cg relative wind velocity component 
along the subscripted axis, ft/sec 


Airplane principal axis system 
Retrerence Cartesian axis system 
Airplane altitude, feet 
Rate of ascent, ft/sec 


Airplane angle of attack relative to principal 
axis system, degrees 


Direction coSines in terms of Euler angles 
(see Table I) 


Airplane yaw angle, degrees 


Lagrange generalized coordinate indicating 
G2EGCCEVON OF Spin radius Vector ,. racduens 


Spin rate, radians/sec 


Aileron, rudder and elevator deflections, 
deqmees 


Airplane principal axis system orientation 
Euler angles with respect to the Xy, Yi. 
Zy system, radians 


Airplane orientation in pitch, roll, and 
yaw with respect to the X,. Yy. Zy system, 
radians 


Pitch, roll and yaw rates, radians/sec 


Incremental angular rotations about the xX 
¥,, Z, system axes, radians 
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ie 


Air density, slugs/ft 


Angular rates of rotation about the air- 
plane principal axes, equivalent to p, q, 
and r, radians/sec 






ie INTRODUCTION 


A. FOREWORD AND ACKNOWLEDGEMENTS 

The intent of this work was to develop an alternative to 
the conventional method of computing spin characteristics 
which thus far has not been uniformly successful in pre- 
meeting airplane spin modes. In this effort, the author 
has tried to approach the problem with the capabilities of 
modern computers in mind and devise a method of solving 
directly for spin modes rather than mathematically tracking 
Peep lane motion through all intervening phases of motion 
preceding the spin. If any success is realized through this 
approach as outlined herein, it will be due in no small 
measure to one tireless efforts of Lieutenant R. L. Champoux, 
USN, whose ccatribution. of months of hard work and extensive 
knowledge of computer programming were instrumental in the 
evolution of this approach. Lieutenant Champoux's separate 
work on the computer aspects of this analysis may be found 
ie Ref. 9. The author is also indebted to Professors L. V. 
Senmidt, T. H. Gawain, and M. H. Redlin, of the U. S. Naval 
Postgraduate School whose counsel helped shape many aspects 


em chnis work. 


B. BACKGROUND 
Airplane spin or autogyration is a mode of motion 
exhibited by virtually all fixed wing aircraft. It is at 


best an undesirable flight maneuver and, with growing 





frequency in new airplane designs, a terminal maneuver from 
which recovery cannot be effected. Although the flight con- 
@ations which lead to Spins can be avoided, pilots of military 
iont attack airplanes and especially fighter aixvplanes do 

not have that option if pressed into flying to the "edge of 
m@@e airplane's flight envelope” by the mission or a tactical 
Ssltuation. 

The need to utilize the edge of the flight envelope in a 
combat role is readily apparent when one considers that it is 
mmey in this fringe area that the Superiority of one airplane- 
pilot combination over another can be realized. However, less 
Obvious is the extent of additional exposure necessary for 
Mmiertning in this area if a pilot is to become proficient in 
meer ticular airplane under such conditions. Thus, flight 
near the edge of the envelope must be more than an occasional 
event for the military pilot. Unfortunately, loss rates due 
tO irrecoverable spin modes in advanced tactical jet airplanes 
meeaecntly makes effective training in this area economically 
Meentbitive. Thus the margin of Superiority built into the 
airplane at great expense is in many cases lost. 

He 1S in recognition of this fact, and the fact that spin 
predictions based on conventional analysis techniques have 
not shown much success, that this effort was undertaken. It 
is not expected that this analysis of airplane Spin will lead 
directly to solutions of the spin problem itself, but it is 
hoped that this work will cast the problem in a new mathe- 


matical perspective, from which additional insight into the 





Spin problem might be gained towards a goal of eventually 


Besigning out irrecoverable Spin modes. 


C. GENERAL 

Historically airplane spin studies have employed forms of 
Ene same equations of motion which are used in the small 
Meeeutbation approach to airplane dynamics. This method of 
mathematically modeling the spin appeared promising because 
moe Equations could be readily integrated by numerical 
methods to produce time-histories of characteristic parameters 
which could then be compared directly with flight test time- 
histories. It was also apparent that the body axis coordinate 
system used in the classical approach provided simpler 
€@uations for random tumbling motion than those obtained 
uSing an inertial system of coordinates. Thus the post- 
Seeul and incipient Spin gyrations of the airplane, which in 
Memy cases are of equal concern, can be more readily tracked 
by a computer in a body axis system than an inertial system. 
moe cegree of success realized by thus extending the classi- 
Cal approach however depends to a great extent on the 
judgment and intuition of the investigator in leaving just 
meee COUpling in the equations to permit Satisfactory 
memeaesentation of the motion without introducing excessive 
mem-linearities. Frequently this process degenerated to a 
“trial and error parametric approach without much insight 
being gained and without general validity being obtained for 


oeatao lane COnmmgquratLons Or Spin mo@es. 
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On the other hand, the treatment of the spin problem in 
an inertial coordinate system appears to add undue complexity 
to the computer tracking problem, and indeed for random 
fumbling motion the non-linearities become excessive. However, 
if the steady, or quasi-steady* Spin is to be studied, an 
inertial cylindrical system offers a number of distinct 
advantages, among which are the ability to differentiate 
between motion due to spinning and motion due to tumbling“, 
@mad a Sufficient reduction in the complexity of the coupling 
Benepermit computer solution of the exact equations>. Since 
the intent of this study was to develop a generalized method 
Of solving for steady spin modes, the equations of motion 
were derived in an inertial coordinate system. The details 
of this coordinate system are shown in Figure 1, and discussed 


in some detail in Section III. 


Spinning motion in which state variables have a mean com- 
ponent and an oscillatory component. 


= Random motion inconsistent with steady or quasi-steady spin 
motlon . 


3No small quantity approximations are made in either the 
demivation Or Simplistiication of the equati onsSaaaaweE on. 


dial 





II. SCOPE 


This work was preceded by a survey of current airplane 
spin literature [Ref. 1] to determine the present state-of- 
the-art in airplane Spin analysis. From that study it was 
evident that possibly another avenue of approach to the 
problem ought to be explored. Thus, rather than attacking 
@ge Spin problem on the broad formidable front of mathe- 
matically following an airplane through all transient 
gyrations from stall to spin, it was decided to limit the 
Seope Of the investigation and probe directly into the final 
Secady state motion of the spin where a force and moment 
Sema librium situation might yield to simpler mathematical 
meresentation. In this approach no restrictions were 
imposed on the form of the airplane aerodynamic data to be 
meeod in this analysis other than requiring that it be ref- 
emenced to the airplane's center of gravity and principal 
axes. By leaving the details of the makeup of such data 
somewhat general, the analysis is freed from any incumber- 
ances of classical aerodynamic derivatives and thus can 
Ppeeceeaymore directly to the derivation of the general 
equations of motion. Although conventional force and moment 
derivatives are suggested for use in initial computer pro- 
gramming, subsequent work may suggest a more appropriate form 
of aerodynamic force and moment data, or possibly more 


representative aerodynamic derivatives may be devised. 
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III. DERIVATION OF THE EQUATIONS OF MOTION 


The coordinate system used in the derivation Cleene 
Sauations of motionseis a combination of an inertial 
Syilindrical system and the airplane principal axis system. 
The vertical (z) axis of the cylindrical system corresponds 
to the central axis of the motion or "Spin axis". The 
advantage of the cylindrical system is that it locates the 
airplane center of gravity in convenient terms of an alti- 
tude coordiante (20), Spin radius coordinate (R) and angular 
position coordinate (%), when the airplane is in a steady- 
State Spin. The orientation of the airplane about its 
center of gravity with respect to the inertial cylindrical 
maecem 1S Specified in terms of Euler angles. A cartesian 
coordinate system is fixed at the cg position on the (R) 
vector having its (x, ) axis in the (-R) direction and its 
(z,) Ss parame! to and un the oppesite direction of (Zo). 
This cartesian system provides a reference for the orientation 
Euler angles. Thus for zero values of the Euler Angles (0), 
(VY), and (d), the airplane would be upright, wings level, 
and have the (X) axis of its orincipal axis system pointed 
inward, along the (R) vector, directly at the spin axis. 

The details of this hybrid coordinate system are shown in 
Fig. 1l along with the positive directions of the position 
coordinates. A general schematic of Euler angle relation- 
mms 15 Contained im Fig. 2 , and Equations for the 


direction cosines of the airplane axes expressed in terms of 
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MmeRTIAL REF. 


Fig. 1. The Coordinate System Used to Model the Spin 


Z»*9O1z 
: ee 


Angular relation- 
Ships between the 
airplane X,y, z and 
the mY Zz, axis 
systems (shown 





Yi 
inverted to 
macilitate ~ 
visualization) X,Y, Plane 
*A negative value 
i as of $ has been 
1 oe Coon depicted in this 
X, 7 ee are schematic 
f 
as (INTERSECTION OF 
X.-Y%. AND x-Y 
X N PLANES) 
Bag. 2. 


A General ecenematic of Euler Angle Relatronsnips 
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Euler angles are listed in table I. The rationale behind 
this choice of coordinate systems iS that 12 mone yeimply 

(in a mathematical sense) represents the motion being 

modeled. Thus even for oscillatory spins the motion would 

be readily recognizable and translatable into mean values 

of the six coordinates (three position and three orientation). 
In addition, Since the principal objective is to solve for 
steady spin modes, the complexity of stall and departure 


motion in this coordinate system has no impact on the results. 


Table I: Direction Cosines in Terms of Euler Aneless 


Cosines of Angles between X, Y,Z and X,, Y,, 2, 








ee ee 
= cos ¢ cos @, = —singcosy 
Xx, : 2 ¥ ai ce sin @ sin y 
: — sing siny cosé@ — cos ¢ siny cos @ 
. & ° s 
= cos ¢ sin @oo = —sing siny ae ek 
Y, i 4 sa @3.9 = —siné cosy 
: + sing cosy cos 6 + cos ¢ cos ¥ cos ¢é = 
a | a3 = sing sing @3 = siné@cos¢ @33 = cosé 


With the coordinate system defined, the equations of 
motion may be expressed in terms of the six generalized 
coordinates z., R, xy , 6, Y and @ by the use of Lagrangian 
mechanics. The only simplifying assumption made in the 
derivation of the equations is that the airplane develops 


meso thrust in the Spin. This is frequently the case in an 


ererence 2, page 158. 






actual spin situation due to jet engine stalls and/or flame- 
outs resulting from inlet airflow distortion. However, the 
restriction imposed by this assumption can be easily removed 
by including negative dissipative work Eexms (On sth osmroat 
hand sides of the Lagrange equations which follow. Such 
Berms would include X, Y and Z force components as well as 
force and moment derivatives about all three axes due to & 
and - xX and 4, and the angular rates p, q, and r. 

pace derivation Of Ehe Gquatlons Of moOGlOon iS sest se iaiuce. 


trated by first writing the general Lagrange equation: 


oT 
£27 ce oT + et =F , v= [+++ 6 (1) 
dro O4r {r tr 
Where: 

T = system kinetic energy 

V -— System potential energy 

qd, - A generalized coordinate, r = 1,2,°°:,6 

Fq. ~ A generalized force due to non-conservative 

LOuCes, Lo" 1,2, °° 40 


For this particular coordinate system, T can be written as: 
T= “aM (ie) 4 (Ry HEP f+ A(T, WO? + Ty? +1202 | (2) 


Where ly, IL. To Oy Wy Wo: are moments of inertia and 
angular rates of rotation about the airplane's principal 
axes. 

Using Figs. 1 and 2, the angular rates relative to the 
mmertial reference can be weteten in terms of the spin 


coordinates as follows: 


16 





G)y = (pry) a3 + 6cos 


Cy = (Wt¥)c23 -OS 


We = (VIE kz + 


(3) 


(4) 


(5) 


mete that # is included to account for the rotation of the 


Ky: Y\° Za system. 


The partial derivatives OT/ 94, for G.=%4,k,0,0,V 


and » SIS e 


OY EMR Yt Tx GX, + Ly ay alz3 + 12Qz 53 


Iy@, cosh -Ly Gy SIN D 


oy x x Xs + Ty Wy k23 + Iz Wz X33 
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(6) 


(7) 


(8) 


(9) 


(10) 


(Ey 





And the time derivatives of these are: 


27) al 
ir (32) = Ma on 


4 (9T) emi (13) 


30p ‘ ' ; 
¢.(S7)= M(2RR¥ +R? Y) 4 Lx (Oo, 3 + Oy X13 ) 
| (14) 
+ Dy (Wy 234 Oy 5) + Lz (rg o(g3 +t Dz X33) 


¢ ($2) - Ix (0), cos ~$ ax SiN ~)-Iy (wy su p+ oy COS $) LS) 


F (BE = 1,(cid 43 + nA) + Ty (dy aes + Gyelrs) 


(16) 
+13 (W, 32 +Ws 33) 

4/2T)- 2 

a($5) = 22% (17) 
Next, the partial derivatives 2 Yq, are obtained: 

oT 

fm =O (18) 
OZo 

er = MYR (19) 
QT . 20 
ay O (20) 
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ae Gx Me | 
ET = 10h SS + yy So + Tez oa (21) 


E 2 Ta, + Ty ody 2% 5 Tz Ode 22% (22) 
Sr = Ta 3 +IyWy 2 a7 +1404 50% (23) 
And from V = MgZ, the values of yt are: 
ov 
= (24) 
366° M9 
OV_dV _OV _wW _D (25) 
OR OF 06° 3p 397° 
Next, evaluating the quantities OO), , OX: 3 and Mis) a, 


eek ~6UcX, UY: Or Zand i ~ We 2, or 3, and r = 4, 5, or©r Ge 


Oks (PY a3 +(Yrt)X3 + Ecos -69SN (26) 
y= (W+¥) 23+ (V4¥) op3 -S HN P- & BcOSP (27) 
= (Y +¥) of33 +(Vt¥)O, +6 (28) 
3 = SCUSOSIND + O5SIN G COS (29) 
X3= 6 COSOCSP-PSING SIND (30) 


de 





X33 = -OSINE 














(31) 
oes Cos 6 SIND (32) 
ee = 31n @ COS> (33) 
Hes = Cos Go Cos (34) 
: = - SIN © SIND (35) 
3 33_ _ SIN © (36) 
dog 0023 _ Doxo3 oD X33 


= = -(0 

op OP TY > TD sci 
The partial ee tee ee etc., are also easily obtained, 

ae 
and will not be shown here. 
Pima, by SUPStituting the derivatives*defimed in 

equations (12) through (25), six Lagrange equations of motion 
Gan be written in terms of the quantities defined in equations 


(26) through (37) and the partial derivatives of the angular 


velocities are as follows: 


MZ, +M¢ = Fe, (38) 
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MR-My'*R= FR (39) 


ME (2RY FRY) 41x (ad 3 + Dy 3) + Ty (Wy X25 + Wy X23) 
+I z (Wyx33+ Wa%53) = Fy _ 


T, (Wy Cos g - $a, SIN ) -Iy (cy sin b+wWyd c0s@ ) 


(Tyla BF + Ty @y OX + Tz@3 Ot) = Fe (41) 


I, (, 3 ft Ws, 3) a Ly (ay 022 + Gy 423) 


ate 2 (a, 33 t Gz X33) = Fy (42) 
2G) 
T, Gz ~( L,.@, 36 + ly cay ) = Fy (43) 


mS diseussed in Ref. 2, the concept of virtual work can 
be employed to derive the non-conservative forces 2 Fe : 
r) 


Fy. Fe. Fy and Fy which constitute the generalized forces 


in the Lagrange eguations (38) through (43). In this 


PomemiLlgque, each of the generalized coordinates (Zo. Rigs 


Q, Y or D ) are increased, one at a time by a small positive 


increment while holding all other coordinates fixed, and the 


expressions for the virtual work terms are determined. 


Since 


system forces and moments due to aerodynamic effects are 


applied at the cg and expressed in orthogonal X, Y coordinates, 


no coupling of the virtual work terms occurs and they can be 


Fadl 





readily derived by inspection using Figs. 1 and 2. fThe 
partial derivative of the sum of the virtual work terms is 
then taken with respect to one of the generalized coordinates 
and the corresponding generalized force is thus obtained. 


For example: 
02, 2 {Fy (gis) ey (gh) aibe, ($2) 
tMy (Oz, (2 ') + My, (9&1) t Mz (3@) fas, Fz, 


(44) 


Where o. , Ey, , and Ex. are positive incremental angular dis- 


placements about the respective axes X ed and Z@ 


cae Ed fe 


eence, bY Inspeceion Of 10g. IL; 


Oz 


fs. (45) 


And, aS a result of orthogonality: 


Ox, _ Oy _o&, _ORy, _ d Ez, * 
Be = 32° 52, * 5a, az, * —_ 


Thus we obtain: 
ae Te, al 


And, by identical procedures: 


Fe = = Sy (48) 
BA Re Mz, aes 
Fe = My Cos y+ My, SIN YW = My, cosd tly Silo (50) 
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Fig Mx, SINE SIND ~ My, SIN 6 Cos Y 
+Mz,Cos @ = Mz > 


Fy = Mz, | (52) 


mae by employing the direction cosine relationships listed 
in Table I these forces and moments can be obtained in terms 


of forces and moments along and about the airplane's principal 


axes. Thus equations (47) through (52) can be rewritten as: 
/z, = - (3 F +23Fy +033 Fe) 2) 
fie (06, Bet Xo) Fy +3, Fa) (54) 


Fe = Ria 2 Fy +%52 Fg) + (Oh My teh ag ty taza) (55) 


if 


Fo = Mx COSD - My SIND (56) 


xe) 





Fy, = (Xs My tes My + %33 Mz ) (58) 


Finally the complete dynamic equations governing motion in the 
chosen coordinate system can be written in terms of the six 
generalized coordinates z., R, 6, 9,Y, and @ by equating 
the left hand terms of equations (38) through (43) with the 
respective right hand terms of eguations (53) through (58). 


These equations are written in their entirety in Appendix A. 
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IV. INCORPORATION OF CONVENTIONAL AERODYNAMIC DATA 


A. DERIVATION OF REFERENCE PARAMETERS 
In order to utilize conventional aerodynamic force and 

moment coefficient data in the equations of motion just 
derived it is necessary to obtain such typical reference 
parameters as airplane angle of attack (®), sideslip angle 
(7) , center of gravity velocity (Veg ), roll rate (P) ; 
meech rate (9), and yaw rate (/) in terms of equation 
variables and their time derivatives. This can be done with 


the aid of Figs. 1 and 3 as follows: 


wt ———_ 


Veg= XE +7, @ + # C3 (59) 


where @,, €,, and C, are unit vectors along the Gy 
25 exes respectively. 


Bue from Big. lL: 


Thus: 
Vege -RO -¥R &z -Z, 3 (60) 
Waerefore the relative wind at the airplane cg is: 


Vew = & & + 0RE, +2, E, (61) 


By wee of the dixection cosines given in Table I, the relative 
wind can be resolved into components along the airplane's 


X, Y and Z axes as follows: 


Vewy = Xj, R+ X)2 ¥R + X23 E, ee 
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Vay = 2, R +X22YR+>X23%, (63) 


View = XR +X, 4 RIAX432, om 






* 
Y 





V V 
tan of = Ruz tan Z = Ruy 
RWX 


Fig. 3 Definition of Positive a and ie! ANnGLes 


Therefore & and /@ can be written as: 


I (hai R thse FR + 033% 
Hy) R 4)2¥R+K)2 Z, f 


of = TaN 


(65) 


atl 1 Mi adhe 
Ky KR F002 R402 2, 


ime magnitude of the relative wind can be written from 


; a he 
B= Tan o2, R +220R M2 ee 


equation (61) as: 


: 7 Re 
Vey = at os (¥R)* + Ae; it 
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Finally »othesrotational raees Pp, Gemeand taltaveut Ene earn. 
plane's x, yy, and Z axes are adentical teste aig acs 
rates do), ; Oy and Co); given by equations (3), (4), and 


(5), and can be written directly as: 


P=(Vt¥ )oy3 + 8 COSD Mee) 
q= (Vt Joes - 51ND (69) 


r= (Vth )oxs3 +p (70) 


B. MODELING THE FORCES FOR THE STEADY SPIN 

| Tabulated aerodynamic data, although virtually always 
appearing in the form of dim@nsionless coefficients, varies 
widely in terms of the reference axis system, point of 
application of forces and moments, and the characteristic 
fmgth used in non-dimensionalization, depending on the 
mmeencded application of the data. By the wse of the princi-— 
pal axis system in the foregoing derivation of the equations 
Guemotien, considerable sSimplication was achieved. MThere- 
Heme it appears that conversion of aerodynamic data to the 
principal axis system and cg as references would generally 
Pemeoceterable to the alternative Of adapting the e€qmacions 
E> the particular form of "Bae data. ~ This "should present ne 


Eemmous problems, however, Simce the equations must be 


Solved numerically and conversion of thewdata by computer 


Zu 






is by comparison a relatively minor additional task. Since 
some data manipulation will likely be reguired in any event, 
the mean aerodynamic chord ree has been used in the £ollow- 
ing development as a common characteristic length in an 
attempt to be consistent and also allow the use of a single 
“characteristic" length parameter as a computer program 
input. However, the equations which follow could be adjusted 
ESD LOpcCicdtely {Ouse Othem echaracreristie lengihs- 

The general form of the aerodynamic forces and monents 
contained in the right hand terms @f equati##ens (53) through 


'o) Can be expressed as) follows: 


Fi = 72 Se. CV7S = ie XY, ore 


Mi = /2 Cu; PV7SC 

Where C has been used aS a characteristic length in non- 
dimensionalizing all moment coefficients. The force and 
moment coefficients ( Cr; and Cu; ) are, in turn, of the 


ZolLLOwing general form: 


Ce; = Ce, + Chin F(R + Gig H4)+ Ge FCP) 
ate Cre. ny, + Crig -F (Sr) Ge : f (Sp) 


Cui =CMi, + Cuz ‘f (p)+ Chai” f(4)+ Cu;4:t(r) 
+ Cuis'F (S«) + Cnig f(r) + Gai gf( Ss) 


Where: 


1. All Crij and CM jj (j))> 1, 2)7°*" 7 andi sey Ore) 


are aerodynamic force and moment derivatives, as 
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defined in Table II below, which must be experi-— 


mentally determined and tabulated for a wide range 
of oC and (3 EUS - 


Table II: Definitions of Aerodynamic Derivatives” 


SUBSERIPT j 





De fp), # (4), and f(r) ate Functions Of oll) peren 
and yaw rate respectively. 


5. & nee and iy are defined as: 
dq = Addlleron deflection 


br - Rudder deflection 


Sh 


A good example of the form of tabular data that could be used 


Elevation deflection 


in this manner iS contained in Ref. 3. 


At this point some comments are appropriate to insure 
that the above method of handling the aerodynamic forces 
and monents are put in proper perspective. First of all it 
Should be noted that the validity of any calculations of 


Spin parameters obtained by this approach to the spin 


Aerodynamic derivatives Wien Geepect to X and 3 have not 


been included Since of and @ are zero for the steady spin. 
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problem depends entirely on the modeling) of the acwegynanic 
forces and monents since the equations of motion as derived 

am Section EEL are exact. See@eondly it shouldme be reagan zed 
that the aerodynamics of the spin problem involve low speed 
three dimensional stalled flow, and as such forces and moments 
are nonlinear. Thus the Superposition of the separate 

Sercects Of trlaxial motatwon, control deflection, wand vether 
Static wind tunnel data cannot be justified analytically. 

The only wind tunnel test technigue that presently appears 

to offer any possibility of yielding satisfactory aerodynamic 
data for spin is the rotary balance force and moment measuring 
apparatus which iS reported on in Ref. 4. 

Next, allowing that for want of better data the conven- 
tional approach outlined above will be used to exercise 
computer program and look for a zeroth order approximation 
femooin Characteristics, there is another problem that stems 
from this application of conventional data. The problem is 
that existing data is generally of insufficient range in &Y 
and 3 to accommodate the high angles of attack and side- 
mero angles characteristic of autogyration. Alse many f 
mien basic set of 42 aerodynamic coefficients suggested by 
the above development describe coupling effects which in 
memvenelonal airplanes at controlled flight altitudes are, 
air general, not Significant and therefore have not been 
measured. Thus, for example, roll moment due to elevator 
aeflection is normally ignored, but may be cf Some Siguata— 


cance ina spinning F4 airplane where the strong negative 
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dihedral of the stabalator could cause the upwind panel to 


act as a rudder with prospin deflection at highoY and 3 


emgles., Certainly in Such Cases more et tChevacrodymane 


derivatives including the 42 above can be ignored without 


testing. 


Now admitting to the fact that force modeling with con- 


ventional stability derivatives will be done with at best 


Mmrneconplete Set of data, it then becomes necessary £o 


exercise caution to insure that those coefficients that are 


available and used will allow a force and moment balance 
be somehow achieved. Thus, if at least one damping term 
not included for every driving term in each of the six 
degrees of freedom, obviously no stabilized solution is 
meeainaole. Finally, in fairness to this effort it must 
pointed out that the modeling of aerodynamic forces and 
memes On Ene Spinning airplane, although recognized as 
the key to the problem under study, was non-the-less 
peripheral to the expressed purpose of this preliminary 


feen. Once a workable computer program has been written 


emamehecked by use of conventional aerodynamic data, sub-— 


sequent efforts should be directed towards improving the 


modeling of aerodynamic forces and moments. 


ok 


CO 


is 
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V. SOLUTION OF sTHE,EOULTUONS (ObSMOTEGN 


A. GENERAL CONCEPT 

From an inspection of the equations of motion contained 
in Appendix A (in abbreviated format) it is evident that, 
although they are ordinary second order differential 
Pemations, they are also so extensively coupled and shel th 
linear that direct analytical integration is out of the 
feisccion. Even if the left hand sides of the equations 
ewe linear, the forces and moments which constitute the 
right hand sides are functions of a three dimensional flow 
mo lvVing the geometry of the airplane. Consequently any 
Semuciton scheme will of necessity involve some form of 
iterative numerical procedure to circumvent direct inte- 
gration. The problem thus becomes one of finding an 
iterative scheme which provides some degree of uniform con- 
Memegence to Steady or quasi-steady values of the six Spin 
parameters ( Z,, K, y a, YW ' and 9 ) without losing 
numerical significance or uSing excessive computer time. 

If the problem iS visualized aS a trial and error 
poSitioning and orientation problem where the airplane is 
ewetcementally reoriented angularly in response to an aero- 
dynamic moment inbalance and Simultaneously adjusted in Z,, 
y and Rin KEeSponse £6 an acroudynam Cc LOnece Inealanec, re 
then appears to be reducible to a simple series of stepwise 
HIyUStMeNtS Or iterations. AlSe, Since the problem woul 


thus be restarted on each successive iteration with updated 
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force and moment data corresponding to the adjusted spin 
parameters, such a procedure would avoid the possibility of 
liaosing numerical significance. The problem Oiadererniinin. 
consistent incremental changes in each spin parameter to 
achieve some degree of uniform convergence of such an 
iterative scheme can be easily handled by uSing the principle 
of conservation of momentum. By rearranging the equations 

of motion in such a way as to facilitate the computation of 
accelerations in each of the Six spin coordinates it is then 
possible to apply a common increment of time and thereby 
compute a dynamically consistent set of incremental changes 
in the parameters and obtain new parameter values for the 
Subsequent iteration. This, ina sense, amounts to follow- 
ing the motion of the airplane as is done in the conventional 
moving coordinate system approach to the problem, but the 
hybrid inertial coordinate system has several distinct 
advantages in this application, as discussed below. 

The rationale upon which this suggested solution techni- 
que is based is that it duplicates in a computational sense, 
what actually happens to an airplane under such flight con- 
ditions. However one key advantage that the computational 
scheme has over the actual airplane is that, on each restart 
Or iteration we can discard any angular momentum inconsistent, 
with a steady spin situation. Thus, we can avoid the problem 
of carrying along tumbling momentum which could considerably 
lengthen the computational time or cause the solution to 


pass right through a steady spin mode or recovery without 
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Stabilizing. Of “course to vbeatiZeuemis advanrage meas 
necessary to be able to differentiate between that part 

of the motion which constitutes tumbling and that part which 
is consistent with steady spinning. Herein lies one of the 
major advantages of writing the equations of motion in an 
inertial cylindrical reference frame rather than in a body- 
eexed moving coordinate system, for in the hybrid inertial 
Secem tumbling is readily recognizable as terms involving 
6, org. 

B. SIMPLIFICATION OF THE EQUATIONS OF MOTION 

Wie ce above approach in mind the equations ©f Mmorlen 
Saminow be Simplified Ssuffieiently to meet the requirements 
Of the solution scheme without the necessity of making any 
restrictive simplifying assumptions. A1so since the 
mieenaded iterative technique deals with the problem directly 
in dynamical units, we do not have to non-dimensionalize as. 
would be necessary if one resorted to a somewhat more 
Bessctract mathematical method of varying the spin parameters 
Mem@otaln convergence. As a result one can avoid obscuring 
m@empnysical nature of the problem and is in a better 
position Lesutilize insight in monitoring, trouble shoecrimea, 
me Lncerpreting iteration trends. 

Mie simplification iS "started by recognizing that only 
steady or quasi-steady spin modes are sought. Thus, by 
setting angular orientation rates ( 6, y ‘ ane) egual to 
zero a considerable number of terms can at once be eliminated. 


The equations contained in Appendix A thus reduce to the 


following form: 
34 








2, Equation: 


ME,+ Mg = — (Cry oa + Coy X23 + Gg X33) OS {(¢K)* 
(73) 
+ (K)* + (%,)* f 
A equation: 
MR-M%7R = 


rs (CE, ny tFy Xz, + Cre 5, ) (re) 


+(K)* +8)? § 


(74) 
¥ Equation: 
MR (ZR Et RE) tT {Uti )0,37 6 Cos Pha eG L(Y 
+9 ys & SING X34 3 {( P+) &,+D} x33 
‘ (75) 
= fE(Guy%3 t Cy 23+ Gia%53)-K(G5%2 


‘ ~ 4 2 
Cy oat egy, f 2 {(R)' + (R)'+ 8) J 


a5 


© Equation: 


I, {( (YY) x, + G CoS] Cosg- LLY +Y)K, , 


Pe ‘2 
- E SIN Df SUD 1,47 C086 SIND c, g , 


~Ly i *C0S 6 Cosh Xz5+ Igi* SING H3 


= (Cy C05 9 Guy $v B) PSE Sc4p) (RI +E.) § 


Boequation: 


IR {(vrI)n,s 7 6 cosh} + Dy fl G+d)X23- GING fx23 
+Iaf(0r0) 433 + foly3 (77) 
7 (Coane * Cay %o3 +p 293) OSE I (fa) +(8)+ +(2) *f 


GEquation: 


I, {(v+¥) X337 o} -I,,9° SIMO COSP K 3 

tly SIN 6 SIND Hr 3 (78) 
= 7, ume Oe 

= (cy,) PEE {(¥0) + (@)+(2.)" 


Although R rH a rare parameter which is not Consistent 
with the steady spin assumption, it effects both the values 
of & and p and also the magnitude of the relative wind, and 
enmererore must be retained in the equations. 

Next, tO put the Cquatirons into a vseable om tor ven 


suggested iteration scheme we simply solve each equation 
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for the corresponding acceleration parameter in terms of the 
other variables. The equations then become: 


x Equation: 


ae So Pee 
Z, = -9-(RM 3 +e, Os +6253) $3 {(7R) 


(79) 
+(R)*+(z,)* 
R Equation: 
R= Y 7R- (C60, + Coy Xe + Gey ~5, ) SF { (eR (80) 


+(R) "+ (z,)*} 
¥ Equation: 
. 2: 2 2 2 z zZ a 
f= Sues, Pe Tyas, + Te 03} [Lets +Dy asa +1255 
-(I,-Iy) 8 SIN 6 SIND COSP -Iz Px33-2MRRY 
rf C( Coy %3 4 Cory X23 + Cu X93) RA Ge0%2 + Gy X29 

‘ , - tf 

+ rgolga)} OS { 0)" E+ (2.)"] | 


& Equation: 


(81) 


a Beant! ay 

6 -{t,costo+ 1 sin 6} [-(1-T(G+¥) SW6 sing cosg 
+(I, SIN" p+ Iy Cos*¢ -I; ) Y* SING COSG+ {any Cos (82) 
Chay sing | CSol(¢R)*+ (BR) (2,)° | 
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W Equation: 


a -| . ¢4 
Y =4I,ag 4 ly Xr ¢ (eee ben x3 +I X23 +14%33)¥ 


-(I,-Iy) 656 SIND Cos $-Isz db c0s 6+ [Cmyx)s (83) 
+ Cmy Xp3 +Ga.%3 | se (yk) + (R)*+ (2) J { 
b Equation: 
b =(F=*)F* swe sin 6 cos b -($+¥) cose = 
< 


+ Cig om { (eR) + (KY (2,)"} 


Before proceeding to develop the details of the pec nenen 
iterative scheme of solving these equations, some comments on 
the steady spin assumption should be made. First of all, 
pealy oscillatory spin modes, as opposed to lightly damped 
steady spin modes that may be very slow in convergence 
Should not, in general, result from the use of mean force 
and moment model test date, the usual form of such data. 
Benoit etae oSerllatory nature of a Spin were attribucable 
to periodic vortex shedding from the fuselage or stalled 
aerodynamic surfaces (von Karman Vortex Street phenomenon), 
the use of mean aerodynamic forces and moments should 
Perectively filter out such oscillation from the computer 
solution. If however the oscillations were strictly a 
moult of the dynamics of the problem, such as a lime 
cycle type of oscillation, then the oscillatory nature of 


the spin should appear in the computations as a failure to 
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Genverge to steady Spin paremeter values. [neti temecse tae 
complete equations would have to be used along with a very 
Mell time increment and the Stability derivatmy——4 cue te of 
and 6 to obtain an accurate description of Ene nceren sn ene 
form of a numerical time history. Although this would involve 
peeing the full equations for the corresponding acceleratien 
termS aS waS done in equations (79) through (84) and would 
also lengthen the computation time it would not be a diffi- 
cult extension and might well provide some valuable insight 
into the stability of such oscillatory modes. This form of 
the Spin equations would be especially valuable in studying 
mae t tect of applying oscillatory recovery control inputs in 
@eeoenance with the spin oscillations to obtain an “autopilot 
om recovery" from an otherwise irrecoverable spin mode. 
feeder Wse of the full equations would be to study the 
Stability of spin solutions obtained by other numerical 
methods. However terms involving ay y Faer= bate p would" re= 
introduce the effects of angular momentum into the computations 
ememwould probably cause excessive tumbling if the full 
~aimeerons Were used to initially solve for Spin modes. here-— 
mee Cither the following computational scheme or some other 
Memerical technique for solving non-linear differential 
equations such as the iterative procedure described on page 
270 of Ref. 7 should be employed to initially solve for spin 


modes. 


C. SUGGESTED COMPUTER SOLUTION SCHEME 
Although the equations are rather complex the suggested 


computer iteration scheme is fairly simple. Basically some 
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starting position orf initial “guess must) Seesur or rede 
terms Of altitude rate ws ices Spin rate, orientation 
(6, WY , and D ), along with control positions b, ; $, and On : 
Angle of attack and angle of Sideslip are then computed frem 
equations (65) and (66) and, in turn of and @ are used to 
obtain values for the various aerodynamic derivatives dis- 
missed in Section MB by interpolating from apprepriatre 
tables. The Six aerodynamic force and moment coefficients, 

CF; and Cy; are then computed as outlined in equations (71) 
mea (72) and, with these values in hand, Z, a Can be 
computed directly from equations (79) and (80). The other 

St We i o7 a a 

accelerations J , on Y) > ane D ,; Can be Gemnputccmey <a 
looping routine encompassing equations (81) through (84). 
This looping will be necessary due to the coupling of the 
respective acceleration terms in each of these equations. 
Once consistent values of Y ; 6 ; y , and d have been 
obtained from the looping iteration, incremental changes in 
each of the seven computation parameters can be obtained 


from the following equations: 


@ 


Zo(n+1) 7 Mommy? won (85) 


Ries) = yt At (86) 


Sone 
Kine = Rin * Rie prttk a Dp 


ee 


Cpadius rate, A, is initially set equal to zero. 
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=} y (88) 
Yn+t) Min) “hati 
: 5 at? 
net) KE re (89) 
= Bs y at (90) 
Pint) 7 Yn) “ee 
7 we 
Pens 7 Q,)+ D ae (91) 


me the tulle Cquatlons were tobe Used, obv feuc haat 
would then be necessary to also compute O ; Y ; D ; G ; 


Y , and D as follows: 


Sine) > Gin)” eat | bees 

nei) = = Gq)" Sinan ‘+6 ate wer 
Vora - Ven) +yat fee) 
Hirst) = Vn) + Ya ytt y at" (95) 
Diys:) = Din) + DAt (96) 
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a 4s wat, - 
aap Ay * PlaeyA® +O |Z se 


The computational sequence outlined above is then repeated 
uSing the revised values of the parameters involved and 
Seneinued until some pre-designated convergence Criteria is 
met, such as all accelerations being less than some small 
=Iomlon Valuc or the rates of change of the various 
mreaneters abe Cachi in turn less than Some Small characecer— 
istic values as individually specified. For example, one 
£OOt per second in e OLS feel per Second ein Z,. one degree 
Ber Second in ie SEC: 

The choice of an incremental time probably will have a 
memy Significant effect on the overall success of this 
Mmeeration scheme. The use of a very small value of Z)\twill 
of course provide the smoothest convergence but at a 
Seqmiticant expense in computer time. A large VG On “ene 
Other hand can easily preclude convergence by drastically 
feeeer ang all parameters from iteration to iteration. To 
illustrate this point and also pin down what is meant by 
large and small At it is necessary to step back a bit and 
consider conceptually what is really being done by uSing an 
meerative Solution scheme. AS pointed out Garlier, the 
equations are highly non-linear and extensively coupled, 
and consequently one haS very little insight as to what 
linearizing approximations would be valid or, for that 


matter, over what range of values can the effects of an 
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individual parameter change be Considered Pincamaae tas 
problem is however overcome in thewysuggested peemation 
Seheme by allowing the dynamics of the problem Go pick vibe 
fagnitude of the inesemental changeswin cach@pamameter. 
Thus if the AG used is of such a magnitude as to constitute 
Pmomall pertubation relative to theycharacteristic timenor 
the problem, the time per revolution, we essentially estab- 
lish appropriately sized small pertubations in each of the 
Spin parameters as well. Thus, for example, a Spin rate of 
40° per second would yield a 9 second turn and possibly a 
mood Starting guess for an appropriately sized At would be 
-0O9 seconds. The validity of this concept will not be 
rigorously or otherwise established here, other than to Say 
Pin in general small pertubation techniques have enjoyed 
seme Success in other engineering applications and this 
meeeats tO be a quite natural extension. 
The problem of picking appropriate starting guesses for 

Z,, R ’ y Oe ome YW , and D unfortunately is not as easily 
handled as deciding on an initialdt. Although some rather 
broad guidelines can be stated such as steep nose down 
pitch, large radius and low spin rate for normal erect 
(upright) spins and shallow nose down pitch, Small radius 
ema high Spin rate for flat spins, the actual searching for 
individual modes must proceed on a trial and error basis 
until an analysis of solutions suggests some better pro- 
cedures. Starting guesses that don't lead to spins will of 
course show up as recoveries to linear stalled flight atti- 


tudes or complete fly-away recoveries. 
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Certainly some inSight can be used to pick Z, aS a 
negative (down) and of an order of magnitude consistent 
with stalled flight velocities (T.A.S.) at the altitude 
where the solution is sought. Likewise, a wings level 
meeatude with the nose down and pointed at the spin axis, 
ora., 9 — 80) Y = ~90° and D = +90°, would be a reason- 
fee Starting Orientation for an erect right Spin. Finally, 

¥ shoula ebvaeusly be Chosen votatleso Ehat 4 tewill ot Gooe 
the airplane out of the ball park under the action of 
Semeritugal force before the aerodynamics can come into 
may. One point Should be made however concerning Searches 
for flat spin modes with small radius values. If conven- 
tional force modeling is used, as discussed in Section IV B, 
then in addition to the general lack of validity of adding 
Meme inear acrodynamic forces and moments, there is the 
@empounding of errors due to the fact that free stream Of , 

V/ and GC vary from point to point over the fuselage and 
wings as a function of orientation when radius is reduced to 
the same order of magnitude as fuselage length or wing span. 
Mies fOr Cxample, if the wings are lévéel and the nose is 
pointed at the Spin axes ina fawey tlatocatkeicudes Ene 
Geometry of the problem will yield a significant difference 
mepoth Sideslip velocity, and consequently (3. from nOSe Jee 
mat! for a small Spin radius. Of course aerodynamic data 
for this type of flow condition cannot be obtained in a 


Seakve Wino tunnel ssetup. 
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D. AUXILIARY ANGULAR RELATIONSHIPS 

The particular set of "Euler Angles" used in the 
derivation of equations (79) through (84) is just one of 
Several variations of Such positioning angles which are 
Generally classified as Euler Angles. The principal 
gavantage gained by the use of Such sets of positioning 
angles is that they describe a Single fixed spatial 
Orientation of a rigid body for given values of Sy Y , 
and On and yield a convenient set of equations for the 
Smee ction cosines: with the order of angular rotation 
uueticit in the equations. However unlike the more con- 
ventional pitch, roll and yaw orientation angles which can 
Memmnecadily visualized once the order is specified, 
Orientation in terms of Euler Angles is difficult to viSualize 
without some sketching, which can be readily appreciated by 
refering to Fig. 2. Thus unless the orientation is close 
EO quadrantal a computer output of orientation in Euler 
Angles may be difficult to interpret. 

iowa llevyiaee thiceSituation, rellationemaps between an 
Sreqered set of orientation angles (1) pitch, (2) roll, and 
(3) yaw, and the Euler angles shown in Fig. 2 were 
ferived. The derivation was facilitated by the sketch of 
the spherical triangle relationships between the two 
pagular sets Shown in Fig. & and involves the wse of 
Napier's analogies from Spherical Trigonometry, Ref. 5. 

The use of Napier's analogies avoids the quadrantal 


ambiguity associated with the solution of right spherical 


45 


triangles by other methods. Thus the equations for con- 
verting the ordered angular rotations e, DQ, and v, to the 


Set of Euler angles used in this paper are: 


an (S*) 
D = by -Tan Sran( Eten ES) }-T aN if mee wey J (98) 
para {| -tu'frn Cerra) 


O= = 2Tan Sin [rent { Ton Ee) Tay )} (100) 
SiN [Tart ae } [teu 


zx Z, 


4 
4 
/ 
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Fig. 4 Schematic of Auxiliary Angular Relationships 
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The equatlons (for vobtatalic S Q. and Vy from Euler angles 


aces 


1 Y OE 
a = D+Tan =a (E) sw (% z =I, a a se (101) 


sin ( 2) Cos (192°) 


y 90-6 yp {0-0 
C Be [Eee Joos (7 ad Tan () sin ( = | ‘iat 


ie 90+ 6 G0+6 
Cos (1248) Cos (toe *) 


featf Tan ieee ] 
b= 40-2Tan SIN [I" 1 sin (4228 
R 


sw Pras! f Te) cs Tran HS 


Cos (129) 


Equations (101) through (103) can be converted for computa- 


ekek) 


tions in radians by substituting WV for 90°. 
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Viz “THE INVERSE PROBLEM 


As discussed in Section ILE B the conventionalmeeiodgat 
modeling the aerodynamic forces and moments on the Spinning 
airplane with eee wind tunnel data lacks mathematical 
midity. Therefore it will be mecessary to develop a 
better method of modeling the aerodynamics of the problem 
if the computational approach to predicting spin modes is to 
be successful. The work done with the rotary balance, Ref. 
4 has been one effort in this area. Such work, however, is 
difficult to pursue without some method of obtaining the 
actual forces and moments as a basis for comparison. 
Pipperently in the past the validity of such aerodynamic 
Giewea has been checked by the indirect method of uSing it to 
solve the partially linearized equation and then comparing 
Semeuted Spin modes with known Spin characteristics. This 
pepeoach may allow some estimate of the overall Success of 
wie Gemobined result of the equations and data but doesn £ 
establish the validity of either. A better approach would 
be to use the exact equations derived in Section V in such 
a manner aS to allow comparison of the 6 computed aero- 
dynamic forces and moments with the corresponding forces 
and moments on the actual Spinning airplane. Fortunately 
this can be done for the steady spin case and thus allow 
Some progress to be made in modeling steady or mean values 


of forces and moments. 
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In order to use the equations of motion derived in 
Section V in this application it is first necessary to 
obtain the values of Z., Koo YP ean D for ene aerial 
airplane or free spinning model. For the case of the model 
in the free-spinning wind tunnel discussed in Ref. 6, | 
determination of the above spin parameters can obviously be 
obtained from tunnel speed, and photographic coverage, and 
thus will not be discussed further. For the actual airplane 
the means of getting these parameters is still fairly 
Simple but in some cases less obvious, and consequently 
Pm be Outlined here. First of all, Z, can be determined 
by the pilot altitude/time observation or photo panel data, 
Ge externally by ground station precision radar or optical 
tracking. Likewise ¥ Can be determined either by the eatee 
om trom ground station photo@raphic coverage. Spin radius 


can be obtained by resolving triaxial cg accelerometer data 


ebeng the principal axes and using the following relationships: 


Get 9y* Je =h0 + (YR) i 


== Fy Va a Gy te -l ee 


P—aorane Orientation in terms of Euler Angles Oe Y , and 
P can be determined from values of pitch, roll, and yaw 


rates and the triaxial accelerations as follows: 


, 106 
P=YSIN@ SIND — 
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oh ¥ SING COSP (107) 





r=yCose (108) 
Thus: 
QD = Tan “A | (109) 


(110) 


@ =Tan’'{ 2 5} 


and from Table IL: 


= a Pe Pe = Ix + Jy Xai tda Xs) 


or solving for {7 


wsc0s' el. agen cna G,. + +9749 | 
] (9, Cc0sg- I, SING) +\9 +(g SUP COSO 49, COSY COSE G5 ney (111) 





tag ai { Bat Ae SIND CaS © +97 C0SB C08 6 ~92 ai 
GCOS - Gy | SIN 


Oey rousiyertnere will be some Joss Of accuracy an Elie 
determination of spin radius and orientation if the Spin is 
Mememexactly Steady. However small fluctuations in the 
mebaxial accelerations and pitch roll and yaw rates Can be 
averaged out and should allow computation of representative 
values of R ; é , 7 pare D for the near steady Spin case. 

Next, G€quations (79) through (84) can be adapted to tic 
problem at hand by setting R and all acceleration parameters 


equal EO 2eLOmp Ine eCqudtions thus Geaduce to, eEhe 1ol2ve7 tna: 


one) 


ce Equation: 
wa Me 
PS f(¥2)* +E) $ 


R Equation: 


ZY “RM : 
Cry ) +Cr, Xp) + Cr X38, 


esfreya(ZF} 


rf Bquation: 


= Cr A} 3 t Cry O23 a Ch X33 


RO Cp, %2 Yr Cry A227 Cea%s2.) 


=C (Cay. %)3 + Cuy X23 +Cy.%32) 
© Equation: 


Ba Ix SiN) Q +I C05" Ia) ¥ “SING COS 
PST (ER) + (4,)°f 
= Cy, COS D ca Cry SIND 


uy BouUat LON : 
Cox %3 7 Cay Xo t Cy 3 i 


d Equation: 


_ 2(t.-ly)¥*sin*e@ SING CSP _ > 


: ae A 
PSE {(¥R)+(4,)* ns 


(112) 


(ia) 


(114) 


(115) 


(116) 


(LL 


Eguations (112) through (117) can also be cast into matrix 


mOrmat LOr Girect Computer Computation as rollows: 


ow 





y Zyxy*RM 


Ay X24 31 || CS { (¢R)*+(2,)} 
X12 A22 A227 Cry = O (118) 
eS!) 
3 X73 X33 |I"Fe C5 (ra)*+(4)* 
CosP - SIND O | ig 
«13X23 %33 [Cy | 
, O | Co a | 
(119) 


— 2( Ix sin “bt Iy C0s*h ~Iz)y¥*sin 8 cose | 





CSE §($R)'+(2,)*} 
= O 
_ 2(Ix-Iy) £* 5106 SH OCOSP 
ese {CR} + (2,)* | 


where the fact that Cie =O , aS expressed by equation (116), 
] 


has been used in equation (114) to show that Cry =O in Bine 
| 


Steady Spin. 


In addition to allowing the computation of the force and 
moment coefficients from steady spin parameters the results 


of this development show that CM * CMy and CH rsligis. /9¥Oile 
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necessarily zero in the steady Spin case, contrary to what 
one might have intuitively assumed.) ingtact , stieowceeag, 
Spin situation where Cu My , and CMe are all zero is a 
special case not likely to be encountered in airplanes of 
Sonventronaledesmen= To aliustracce Ents (point, cauaruens 
(115) and (117) can be solved for the case where the moment 


coefficients are all zero with the result that: 
z Z 
I, SIN'*Y+ Ty CosH=I, (120) 


and 


L, =Ly (121) 


or Simply: 
I, = Ly =Iz (22) 


Meus 1t 1S seen that in a conventional airplane where it is 
unrealistic to have equation (122) satisfied, aerodynamic 
moments must be present to balance the gyroscopic moments 
and gyroscopically precess the inertial orientation of the 
aurplane sufficiently to maintain a fixed orientation 


relative to the spin axis and radius vector. 
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VIE. RESULTS 2D DISCUSSION 


@ ANALYTICAL, ASPECTS 

PAthouGgh Ene pmalfary intent of Ehisseffortewas to 
develop another method of computing airplane spin modes 
from wind tunnel test data, the equations of motion derived 
mee section IV can also be used to advantage in a qualitative 
mm vecal approach to the problem of determining Spin 
SmWaracteristics. In addition to possibly suggesting new 
Peeametric indices of spin characteristics, these equations 
Should also allow a number of experimentally observed effects 
Semparameter variations to be put on a firmer analytical 
foundation, such as the effect of the sign of Bh on 
Spin recovery [Ref. 8]. The fact that the inertia character- 
istics appear so prominently in equations (79) through (84) 
confirms test results which have suggested that the gyro- 
Beopme: aSpeets off thesproblem are as important as «he 
mame vhantes. Another result of this analysis o£ the Spam 
PpuookLem Was to emphaSmze the non-linearity of the sacro- 
Geynamics associated with the spin. Thus computational 
Eeameniques Which use Superposition in medelang the aero- 
dynamics or employ partially linearized equations of motion 
Should be expected to yield rather poor results where high 
solmeeates Or small Spin radii ane invobyed. These are 
probably the main reasons why the conventional approach has 


not really been successful in treating spins in general, 
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and has been particularly unproductive in predicting lar 


Spin, modes: 


B. THE COMPUTER PROGRAM 

Although the suggested computer program outlined in 
Section V C uses conventional aerodynamic modeling which is 
meeent rely nmepresentative in the Spin, it is felt that 
development of a computer program capable of solving the 
—aiations Of motion 1S the next essential~step in this 
eeeroach to the Spin problem. Once the computer program 
evorhing, improvements in the modeling of the forces and 
mements can be incorporated as they are devised. Thus the 
Computational and aerodynamic parts of the problem can be 
Puecued, in the interim, aS Separate problems. It is with 
this in mind that the author has chosen to discuss a few of 
the aspects of the computer solution problem rather than 
address the aerodynamics of the problem. 

Miecwsitcerative scheme outlined in Section VB may exhibae 
a tendency toward a lightly damped oscillatory form of con- 
Pergence Since it is essentially following the dynamical 
Meron Of the airplane, which is characteristically lightly 
Gemmeed in a Spin. The oscillatory nature of the convergence 
is compounded somewhat because we cannot be sure how far a 
Starting guess may be from a force and moment equilibrium. 
To improve this situation there are a number of ways to 
provide some synthetic damping to the iteration to speed 


convergence. Some of these ways are: 


a. LE excursions in radius preve tor be a Deo bleinea 
a result of momentum in radius rate (R), the motion can 


be damped as follows; 


Rene = 4 Kin) + RAT (23) 


where 4 is some factor between zero and one. 
b. if Excursions in altitude rate become a preblempracy 


can be handled as follows; 
e -~ o / ce m 
2 (nu) = B(n) BE AC — 


POO on Pooler Ssthat nay sandse. 1S Cha G Convercenee milan 
orientation may be deadbeat, depending on the size of AT. 
fms 1S due to the fact that we have eliminated jorientation 
meeaees to reduce any tumbling tendencies. Thus it may be 
feeessany tO augment the rage of change in orientation eer 
Meematlion to speed convergence. This can be done in the 


ZeLlLowing manner: 








Bat (125) 
ones) On) + b ( at } 

VATE 
Yas = Yay * © (Fe 3 


| Bart 
Pint) = Yn) + b( Z ) — 


Where bis some factor greater than one: 
Paiciet) aoe col => Gace ic a leas in ¥ may be out of the 


ball park and cause large changes in radius prior to reaching 


a6 





a magnitude consistent with the aerodynamics of the problem, 
it may be necessary to augment the rate of change of this 
characteristically unresponsive parameter. This may be 


accomplished in the following way: 


é 


Yn) = Yn) + C ( yt) (128) 


where C is some factor greater than one. 


©. FOLLOW-ON STUDIES 

As mentioned earlier, the equations of motion as derived 
mm this paper should make it possible to provide an analyti- 
cal basis for many of the heretofore experimentally observed 
Smatracteristics of airplane Spinning motion. Such a cor- 
relation study would be relatively inexpensive to pursue and 
would probably provide a wealth of additional insight into 
the spin problem as a whole and possibly be beneficial to 
the more immediate goals of computer program development and 
meree-moment modeling. It might also be instructive to 
erogram the full equations of motion, as contained in the 
aependax, on an analog computer having a digital computer 
interface capability for the aerodynamic data. Some other 
aspects of the spin problem which could be pursued once the 
computer programming and force modeling problems are solved 
are: | 

(1) Determine the spin modes of currently operated 
tactical military airplanes and investigate the effective- 


Ness OF Lecommemeacd Spin recovery techniques. 


5) 





(2) Determine what characteristic aerodynamic parameters 
can be most economically fixed in existing airplanes with 
irrecoverable spins and investigate the magnitudes of such 
changes that would be required to make recovery possible. 

An example would be providing rudder stop over-rides or 
@aditional deflection for Spin recoveries. 

(3) Study the feasibility of providing strap-on devices 
such aS a deployable fabric ventral fins or wing mounted 
metro-rocket pods for use on tactics training flights. The 
MmEeen Ot Ene adevice Or necessary thrust requirements could 
be fairly easily determined with the assistance of a working 
computer program, or, in the case of the retro-rocket, could 
also be determined by the "Inverse Problem" Technique out- 
immed in Section VI. 

(4) Determine “design envelopes" that prescribe limits 
fiijecerms Of common aerodynamic derivatives in order to 
Pesist the designer in avoiding troublesome combinations 
wmenrever feasible. Such work would also provide the basis 
Pm 1MpLovementsS in the current Military Specification 
MIL F-8785 (ASG). 

(5) Study the feasibility of utilizing an autopilot 
meradpply oscillatory control inputs in resonance with the 
Gynamics of the spin to break an otherwise irrecoverable 
mode. 

These are only a few applications that might be under- 
taken to alleviate the spin problems inherent in the 


designs of todays high performance tactical military 
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airplanes. However the satisfactory accomplishment of just 


a few is all that is really needed. 


a 


VIII. CONCLUSIONS AND RECOMMENDATIONS 


The approach to the Span problem outlined in this paper 
Should provide a means of obtaining some additional insight 
mee the problemeot alrcraft spin. Once a successful 
Gomputer program hes been written and a more representative 
means of modeling aerodynamic forces and moments has been 
@evised it Should be possible to Solve the equations of 
motion and obtain some idea of the spin modes that a given 
airplane deSign may exhibit. The results of Section VI 
mircrcate that gqyroScopic moments play an important rolkevin 
the spin problem. Thus it may also be worthwhile to treat 
the spinning airplane as a Slow turning free gyroscope, of 
betatively high inertia, being acted upon by triaxial aecro- 
dynamic forces and moments. Such a formulation of the Sor 
would allow the application of the broad scope of knowl- 
eeee Of gyroscopic motion from the fields of physics and 
Mechanics to the airplane Spin problem. In thiS manner it 
may be possible to obtain some additional understanding of 
M@emstability Of various Spin modes and possibly avoid some 
Mincesibtaole airplane Spin characteristics without having 


Getailed knowledge of the aerodynamics involved. 
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APPENDIX 


SPIN EQUATIONS 
The foitlowing equations constituce the full ecquationcwom 


motion for a spinning airplane as written in a cylindrical 


Soordinate SyStem. 


Z Equation 


at ead Z 5 wh 
MZ, +MQ = = (Cp 134 Oey X03 + eX) ES 4 (78) +(e) +(Z,) 


K Equation 


LS Ca, we (ee 
MR-MNY “R= -(&, Oy tCey Xo + Cr, ) J {(¢aP+(®) *(@,) 
X Equation 


me (2b ¥+R¥) +Ly { ( Yr Fa .+(p+H)(é CoS @ SIND 
+DSING Cosh) + 6 COS -SPSIN GD f X13 4 Tx I (y+¥)oe,5 
+6 cos flecsesind+dsno cos g f+ Td Wi+¥) og 3 
+(U+¥)(6 ws Cos - PsN sw) -S srg 

EP cos b fevgg +Ty [(V+d)wz3 - Es } (6 C080 caSp 
-g5n6 sf +1z E (G+ ¥)ov23+ (par)(-6 sino) 

+B fotsg + Te {(Wrt)oy3 +O} {-6 smo 
= £E (Luong + Cy 23+ Cua 133) R(Gex %2 + Coy 822 

+ Ces Xz2) | cS {(/0™ RI) f 





S Equation 


LA (U+F aj; + WT) 6C056 sh + sW@ cosh) 

+ S cash -6¢6 smb} CSO -I1, b [GL ) 043 +6 C05} np 
-ly {(U4¥ Jokes +(U+V)(6 0056 Cosd-$5INO SIND) 

-E sind ~6P Cds ofS b-1,¢ { (W+¥) 23 

-ESING feos as {(b+¥)o%3 +6 CoSd} A (p 

+¥ )C0S 6 SIN p} -Iy {(p+y)%23 -6 sing} f(y 

+¥) cose cosh} + Ie {Ut¥) asa +p} {(prt)sm 6} 

= (Cyy CoS P ~ Cuy So) SOLE (a4 (é,)*} 


Y Equation 


Z, { (GY oly + ie cose suip+ sing caso )+ 8 cas 
-E6d wn G fx; +1 L.A (O43 +6 wsg} {6 CSE SING 
+P SINE as, {+ F oo, + (U+8)(6 6056 cos 
-B SIN 6 SIND)-6 SING - rr me: {(g 
4) 0,3 - GSN D}LS (OSG COSP-PSNG SING} 

tLe {(Ot¥ )ly3 + (t8)(-6 NG) + 9 fg3 + {(V 
+) 33 +bff-6 5m 6} 

= (Ooty 3 + Gayl tue tay ) SE fix)? + (P+ (27° 
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GP Equation 

Z, f (G4 ¥) «33 + (Ute \(-Esw 6) + p } Lx {(p+¥ ex) 
+6005 bf] (UY) Su 6 C05b- 6.500} -Iy {G48 aes 
-63N B }4-(UH¥) 50) © Su) gd - Ecos p} 

Gu ) 7g { Cig)? + (8)* +(2,)"} 
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